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Введение. Кандидатская диссертация проф. А.В. Белоконя 
(1941-2013), который на протяжении многих лет был ректо- 
ром Ростовского государственного университета, а затем Пре- 
зидентом Южного федерального университета, была посвя- 
щена асимптотическим методам решения контактных задач 
теории упругости для тел цилиндрической формы (1969 год). 
В настоящей работе исследуется контактная задача теории 
упругости о кручении жестким вкладышем бесконечного по- 
лого цилиндра, внешняя поверхность которого жестко за- 
креплена. Вкладыш находится внутри цилиндра и имеет ко- 
нечную длину. Эта задача была в 1971 году поставлена и изу- 
чена А.В. Белоконем, который при помощи интегрального 
преобразования Фурье свел ее к интегральному уравнению 
относительно неизвестного контактного — напряжения. 
А. В. Белоконь получил полное решение данной задачи для 
случая толстостенных цилиндров, когда символ ядра инте- 
грального уравнения может быть аппроксимирован функцией, 
соответствующей кручению пространства с цилиндрической 
шахтой. В настоящей работе рассматривается преимуще- 
ственно случай тонкостенных цилиндров, дополняющий ре- 
шение А. В. Белоконя. Материалы и методы. Материал ци- 
линдра считается линейно-упругим. При решении задачи ис- 
пользуется метод интегральных преобразований. Для реше- 
ния интегрального уравнения применяется сингулярный 
асимптотический метод. 

Результаты исследования. На основе исследования свойств 
функции-символа ядра интегрального уравнения предложена 
новая специальная легко факторизуемая аппроксимация этого 
символа, пригодная для цилиндра любой толщины. Для 
нахождения оптимальных параметров этой аппроксимации 
используется метод Монте-Карло. Расчеты сделаны в основ- 
ном для тонкостенных цилиндров. В результате в аналитиче- 
ском виде получено асимптотическое решение интегрального 
уравнения. 

Обсуждение и заключение. Найденное новое решение эффек- 
тивно для относительно длинных жестких вкладышей, длина 
которых превышает внутренний диаметр цилиндра. Метод, 
основанный на новой аппроксимации, остается применимым 
и для случаев, когда цилиндр можно рассматривать как ци- 
линдрическую оболочку. Асимптотическое решение можно 
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Введение. Контактные задачи о взаимодействии жесткого бандажа с бесконечным сплошным упругим ци- 
линдром и о нормальном контакте жесткого вкладыша с поверхностью цилиндрической шахты в упругом простран- 
стве рассматривались в работах [1-4]. При этом в сингулярном асимптотическом методе использовалась аппроксима- 
ция символа ядра интегрального уравнения в виде суммы двух легко факторизуемых функций, что позволяет найти 
только приближенное решение уравнения Винера-Хопфа [2]. Позже для этих же задач была предложена аппроксима- 
ция символа ядра одной легко факторизуемой функцией [5,6], что позволило существенно облегчить применение син- 
гулярного асимптотического метода, который получил широкое распространение при решении задач со смешанными 
граничными условиями [7,8]. Анализируя научное наследие профессора А. В. Белоконя [9], сталкиваешься с еще од- 
ной интересной задачей, рассмотренной им в начале 70-х годов прошлого века. Это задача кручения жестким вкла- 
дышем полого цилиндра, которую тогда удалось эффективно решить лишь для толстостенного цилиндра [3, с. 230]. 
Цель настоящего исследования — получить решение этой задачи для тонкостенного цилиндра на основе сингулярно- 
го асимптотического метода и новой эффективной аппроксимации символа ядра интегрального уравнения. 


Материалы и методы. В цилиндрических координатах рассмотрим контактную задачу кручения полого ли- 
нейно-упругого цилиндра, внутренний радиус которого равен К, а внешний А!. Внешняя поверхность цилиндра жест- 
ко закреплена. Внутрь цилиндра помещен жесткий цилиндрический вкладыш длины 2а, к которому приложен крутя- 
щий момент М. Под действием этого момента упругий материал цилиндра в области контакта испытывает угловое 
перемещение и,=-=. Материал цилиндра характеризуется модулем сдвига С. При заданных величинах К, К\, а, = тре- 
буется определить контактное напряжение р(<) в области контакта. Затем может быть определен момент М. При по- 


мощи метода интегрального преобразования Фурье задача сводится к следующему интегральному уравнению ( | <а ) 


[3]: 
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Здесь Г[.(и), К„(и) — модифицированные функции Бесселя [10]. Безразмерный параметр © характеризует тол- 
щину стенок цилиндра. При ©—>›00 функция Г(и) вида (2) стремится к функции 
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соответствующей случаю кручения пространства с цилиндрической шахтой. 
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Ранее было установлено, что при 9>3 и и>0 функция С,(и) вида (3) достаточно хорошо аппроксимирует 


функцию Г[(и) и было получено полное решение уравнения (1) с символом ядра (3) [3]. Отметим, что наибольшее от- 
личие этих функций наблюдается в нуле, где 
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В бесконечности функция (2) имеет асимптотику: 
Ци) - — = -} и >> =. (5) 

и 2и и 
Введем безразмерные обозначения 
К 

в Ио ве (6) 
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Штрихи далее будем опускать. Параметр /, характеризует относительную ширину области контакта. В обо- 
значениях (6) уравнение (1) примет вид 


1 х-< 
ГРОК| == |4 = пв. (7) 
и ^ 
Для решения уравнения (7) применим сингулярный асимптотический метод [5,6], эффективный при доста- 
точно малых значениях /^. 


Результаты исследования. Получение решения задачи основано на сведении уравнения (7) к интегральному 
уравнению Винера-Хопфа, при решении которого используем новую аппроксимацию 
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В 3 
Са оч Е) =А, Е=10. (9) 


Аппроксимация (8), (9) учитывает поведение Г(и) в нуле и бесконечности, см. формулы (4), (5). Кроме того, 
функция (8) легко факторизуема. В табл. 1 для разных © даны значения параметров аппроксимации (8) и ее относи- 
тельная погрешность 0 (в процентах) при и>0. Значения параметра С определялись из условия (9). Значение ©=0 со- 
ответствует аппроксимации предельной функции (3). 

При проведении расчетов параметров аппроксимации использовался метод Монте-Карло. 
Таблица 1 
ТаЫе 1 
Параметры аппроксимации 
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о, Ао В С 0 ©, Ао В С 0 
со 1,970 0,8872 | 0,8447 2,0 1,07 0,4068 | 4,001 8,696 1,07 
3 0,4649 5,651 4,674 0,5 1,06 0,1582 | 6,088 10,25 
2 0,8586 1,150 3,340 0,5 1,05 1,251 5,091 4,580 
1,5 2,192 4,183 3,631 0,3 1,04 1,448 5,002 4,924 
1,1 3,165 6,887 8,878 2 1,03 1,400 6,991 6,876 4,5 
1,09 1,269 6,582 10,09 3 1,02 2,082 5,909 7,875 5 
1,08 2,037 6,583 9,164 3,5 1,01 2,715 6,507 9,500 5,5 








Значение Е специально выбрано большим (см. формулу (9)), поскольку в процессе получения сингулярного 
асимптотического решения используется еще одна аппроксимация в полуплоскости Кез> 0: 


ехр(у(5)) =1+ у(5) ‚ (10) 





= 3 п - Е 
у 21/5? - Е” Е 
При выбранном значении Ё относительная погрешность аппроксимации (10) при действительных 5>0 не 
превышает 0,3 %. 
В результате применения метода Винера-Хопфа главный член асимптотического решения интегрального 
уравнения (7) при малых /^. можно построить в форме 


5 1+х 1-х 1 
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Здесь ег >) — интеграл вероятностей, На основании формулы (11) можно получить выражение для инте- 
гральной характеристики решения (момента М). 
Как показывают расчеты, погрешность асимптотического решения (11) при ^<! не превышает (6+0)%, где 0 
— погрешность аппроксимации (8). 


Обсуждение и заключение. Аппроксимация вида (8) пригодна не только для тонкостенного цилиндра для 
цилиндра любой толщины (при любых значениях параметра 9). Как правило, с уменьшением толщины цилиндра (с 
уменьшением ©) погрешность 6 возрастает. При значениях 9<1,02 цилиндр можно рассматривать не только в рамках 
теории упругости, но и в рамках теории цилиндрических оболочек (толщина оболочки обычно составляет не более 
двух процентов радиуса срединной поверхности). Приближенное аналитическое решение (11) можно рекомендовать 
инженерам для анализа прочностных характеристик деталей цилиндрической формы при кручении жесткой шайбой 
(вкладышем), когда длина шайбы 2а превышает внутренний диаметр цилиндра 2К (при этом А=А/а<1). 
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